Frontier between separability and quantum entanglement in a many spin system 
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We discuss the critical point x c separating the quantum entangled and separable states in two 
series of N spins S in the simple mixed state characterized by the matrix operator p = x\4>X<j)\ + 
jprlni* where x G [0, 1], D — 2S + 1, I d n is the D N x D N unity matrix and |0 > is a special 
entangled state. The cases x = and x — 1 correspond respectively to fully random spins and to 
a fully entangled state. In the first of these series we consider special states \<j> > invariant under 
charge conjugation, that generalizes the N — 2 spin S = \ Einstein- Podolsky-Rosen state, and in the 
second one we consider generalizations of the Weber density matrices. The evaluation of the critical 
point x c was done through bounds coming from the partial transposition method of Peres and the 
conditional nonextensive entropy criterion. Our results suggest the conjecture that whenever the 
bounds coming from both methods coincide the result of x c is the exact one. The results we present 
are relevant for the discussion of quantum computing, teleportation and cryptography. 



Quantum entanglement is a quite amazing physical 
phenomenon. It was first discussed in depth in 1935 
by Einstein, Podolsky and Rosen (EPR) [Q and by 
Schroedinger ||, and since then by many others [p|-p0[. 
This phenomenon has attracted intensive interest in re- 
cent years due to its applications in quantum computa- 
tion, teleportation and cryptography, as well as to its 
connections to quantum chaos. Special importance is 
attributed to the frontier which, in mixed states, sepa- 
rates the quantum entangled from the so called separable 
states. Separability, precisely defined in what follows, is 
essentially a classical concept, and it is the one on which 
the possibility of a description in terms of local realism 
with hidden variables is currently built on. The state of 
a system composed of subsystems A and B is said sep- 
arable (or unentangled) if and only if its corresponding 
density matrix pa+b can be written as 

w w 
Pa+b = Pa ® Pb (Pi > Vz; ^> = 1). (1) 

t=l 1 = 1 

The frontier between quantum entangled and separable 
states for a variety of mixed states is currently discussed 
in the literature. Several necessary conditions for separa- 
bility are already available; however, no general sufficient 
conditions are yet known. Among the existing necessary 
conditions two criteria have shown some preference in 
the literature due to their simplicity and the fact that in 
some cases these conditions coincide with the exact ones. 
Those are the Peres criterium |l3| based on the positivity 
of the partially transposed density matrix, and the gen- 
eralization of standard entropic methods, given by the 
conditional nonextensive entropy criterium, recently ad- 
vanced by Abe and Rajagopal In the present paper 
we are going to apply and compare both methods by ap- 



plying them to several distinct density matrices coming 
from some many body systems. Both criteria recover ex- 
act separable-entangled frontiers in some simple mixed 
states of a two spin 1/2 system |2^j2^], but only pro- 
vide a bound in more complex cases. Although no rig- 
orous proof is available, as we shall see, the analysis of 
some particular cases suggests that generically the partial 
transpose method overestimates the separable region ei- 
ther equally or less than the conditional entropy method. 
Moreover, our results indicate that whenever both meth- 
ods yield the same result, this result is the exact one. 
It is clear that the discussion of quantum entanglement 
in many body systems is crucial both conceptually (e.g., 
for understanding the interface between the classical and 
quantum descriptions of the real world) and practically 
(e.g., for efficiently implementing the input and output 
of a device such as a quantum computer). Some spe- 
cific A-body systems have already been discussed in the 
literature (see ]23|,|2~i| for some classes of systems with ar- 
bitrary A, and [|25[ for the Greenberger-Horne-Zeilinger 
state of N = 3 spins 1/2; sec |2(| for a review). 

Let us consider a A-body system composed by spin-S 1 
operators acting on a Hilbert space H = Hi (8 • • • <8 H n 
of dimension D N = (2S + 1) . An arbitrary density 
matrix acting on this space should have all its eigenval- 
ues non negative. Consider an arbitrary orthogonal ba- 
sis |si, . . . , sn > (si — —S, . . . , S) spanning the Hilbert 
space. The Peres criterion appropriate to this A-body 
systems asserts that from a given density matrix p, with 
elements < sj., . . . , sn\p\s[, . . . , s' N >, it is possible to de- 
fine other density matrices ( j' ll '"' ,ln , where the spin coor- 
dinates i\, . . . , i n . . . , i n — 1, . . • , A), are transposed 
("in <-> out "), i. e., < si • • • s N \a ll '---' in \s[ ■ ■ ■ s' N >= 

< s\ ■ ■ ■ s^-is'^s^+i ■ ■ ■ s- 2 • • • s' in s in+ i ■ ■ ■ S N \ 
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The non-negativity of these new density matrices give 
us the Peres criterion, namely all eigenvalues of a 11 <■■■>**> 
should be non-negative. 

In order to probe and compare the above mentioned 
methods we will consider several sets of density matrices 
of general type 



p(x) = (1 — x)p(0) + xp 



(2) 



where p(0) = px ® . . . pm — 1 d n / D N , 1 d n being the 
D N x D N unity matrix, is a fully separable state and 
p = \<j) X <f)\ \s a. density matrix corresponding to a 
known entangled state \4> >. Then as we increase the 
value of x, there exists a critical value x = x c where 
the quantum entanglement takes place, and we want to 
obtain bounds for x c . 

The matrices er 11 '---' 1 ™ obtained by the partial transpo- 
sition at indexes . . . , i n ), are given by 



(x) 



1 



N 



where 



< 



si,. 



D 



, Sff\a' 



D N 



s li ■ ■ • ! s i t s i 2 ' ' ' i 



. .,sjv|p|si 



(3) 

>=< 

■ S i n 

s' N >. The D N eigenvalues of <r l1 '— ' %n (x) are given by 
Ei = (l-x)/D N + xE l , where E t (i = l,.. .,D L ) are the 
eigenvalues of cr 11 '"' ,Zn . The condition for non-negative 
eigenvalues imply that quantum entanglement is present 
at least for x > x^ (P stands for Peres), where 



(4) 



1 - D N E 



and E is the lowest eigenvalue of cr ll '" -ln . In principle we 
should consider all the possible transpositions i%,...,i n 
but, in all the applications we considered, the transposi- 
tion of a single spin, i. e., tr n , gave us the best (i.e., the 
lowest) bound for x c . 

The second criterion we analyze was introduced re- 
cently by Abe and Rajagopal |2^| . Exploring the classical 
composition of probabilities arising from the nonexten- 
sive entropic form on which nonextensive statistical me- 
chanics {27| is based, they provide a condition for quan- 
tum unentanglement. Within this entropic approach the 
condition for quantum unentanglement is that the con- 
ditional probability for any subsystem forming the whole 
system should be non-negative. In this entropic approach 
the entropy of a iV-spin S system A\ + ■ ■ ■ An, formed 
by subsystems Ai (i = 1, . . . , N), with density pa 1+ ...a n 
is given by 



S q (A x + ---A N ) 



1 



9-1 



(5) 



where q is an arbitrary real number. The particular case 
q — > 1 recovers the standard von Neumann entropy. The 



conditional entropies, taking into account a partial infor- 
mation of the system, are given by 



5,(^i + •• • + At-i + A h+l + ■■■+ A in _i + A ln+1 
■■■ + A N \A n +A l2 + --- + A ln ) = 
SgjAx + ■■■ + An) - S q {A h + A i2 + ■ ■ ■ + A in ) 
1 + (1 - q)S q {A il + A i2 + • •• + A in ) 

where S g (Ai + ■ ■ ■ + A n ) is given by (||) with 



(G) 



^A 1 ,...,A il - U A il+1 ,...,A in _ 1 ,A in+1 ,...,A N [pA 1 +---+A N }- ( 7 ) 

The condition for non-negative conditional probability 
imply that, for arbitrary combinations {ii, . . . , (|^) 
should be non-negative, and consequently the bound for 
x c is evaluated by imposing 



^Pa 1 + -A n = r ftp q Ail +-+A i 



(8) 



The left side of the above equation is simple to calculate 
for the special set of density matrices (||) by going to the 
basis where p is diagonal, and we obtain 



l-x,„ ,l + (D N -l)x 



From (|) and fj|) we obtain 
1 -x 



PA n +---+A, 

where 



D 



N 



D N - n I Dn +x~p AH+ ... +Ai , 



(9) 



(10) 



P~A il ^ VA, n — T^A H ,...,A H -i,A H+ i,...,A in -i,A in+ i...,A N 

[pA 1 +-+A N ], (11) 

is a D n x D n matrix calculated from (||). If we denote 
by Vj (i — 1, . . . , D n ) the eigenvalues of p the relations 
©-© give us 



. 1 — x N 



(^-!)(^lf) 9 + ( 

D n 



1 + (D N — l)x 



D 



N 



i=l 



D 



N 



D N - n +xv % ) q 



(12) 



This equation will give us x c (q) for each value of q. Due 
to the monotonic properties of the entropy (^|), the lowest 
bound for x c is obtained in the limiting case q — > oo. In 
this limit we can neglect the first term in the left hand 
side of (^|) and keep in the sum only the term corre- 
sponding to the largest eigenvalue v(N) of the matrix 
p Ati+ ... +Ain , which gives 



x s = 



(13) 
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where the superscript refers to entropy. In principle we 
should consider the conditional probabilities of arbitrary 
subsystems. In our applications the largest eigenvalue, 
i.e., the one which produces the most restrictive bound, 
is obtained when we consider, in these conditional prob- 
abilities, the subsystem A\ + . . . + Ajv-i< 

We are going to consider initially density matrices 
written in vector basis invariant under charge conjuga- 
tion, or spin-reversal symmetry. In the S^-basis the 
basis vectors are given by the non-null combinations 
\<t>1 >= ■^{\si,...,s N > +c\ - si, . . . , -s N >), (i = 

1, . . . , D N , c = ±1). For S = \ and N = 2 this basis re- 
covers the standard Bell basis. We now consider the den- 
sity matrix p(x) given by (g) where p = \<f>^ >< 4>^\ is 
formed by an arbitrary basis vector with charge conjuga- 



tion eigenvalue c 



■1, i. e, \<j>l >= (|1 > -|2 >)/V2, 
with |1 >= |si, . . . , sjv > and |2 >= | — si, . . . , > 
({si} arbitrary). It is simple to convince ourselves that 
in the standard basis the only non-zero elements of p in 
(f| are given by < l|p|l >=< 2|p|2 >= - < l|p|2 >= 
— < 2|p|l >= 1/2. In order to apply the Peres cri- 
terion we consider the density matrix a defined in (|^) 
by the transposition of the last spin. The only non-zero 
terms of the related matrix a are given by < l\a\l >=< 
2|5j2 >= 1/2, < 3|a|4 >=< l|p|2 >=< 3|cr|4 >=< 
2|p|l >= —1/2, where |3 >= |si, . . . , Sjv-i, — sjv > an d 
|4 >= | — si, . . . , — Sjv_i, sn >• This last matrix has the 
lowest eigenvalue —1/2; consequently, from (Q), we get 
the bound 



x p = 



1 



l + D N /2 



(14) 



for all values of D and N. 

In the application of the entropic bound, the most re- 
strictive condition happens when we consider the condi- 
tional probability of the subsystem composed by (N — 1) 
spins and it is not difficult to see that p in ( |Io|) has only 
two nonzero elements ( = 1/2) in the diagonal, that gives 
the eigenvalue v = 1/2 and from (|l^) the bound 



1 



(15) 



2 D-l 



1+2 



1/2 both 
1 but 



Comparing (|l4j) and (]l5|) we see that for S 
criteria give us the same bound x P — xf = 
for general values of the spin 5^1/2 the bound com- 
ing from Peres criterion is more restrictive than that of 



the entropic criterion since x% < 



It is interesting 



to remark that in the case S = 1/2, where the bounds 
coincide the value is known to be the exact value for x c , 
for N = 2 pi or N > 2 @. The bounds (Juj) give us an 
interesting result. Suppose we need, as would occur in a 
quantum computer, to couple the special entangled state 
\4>k > with a white spectrum environment. The mixed 
system will be described by the density matrix given by 
(^, where the parameter x controls the coupling with 
the environment. The bound (jlj) tell us that we keep 



the whole system entangled as long as the ratio between 
the component of our special vector \<f> k > and any other 
vector in the environment, forming the density matrix, is 
larger than (1 + (D N - l)xf )/(l - xf ) that as N -> oo 
tend towards 3, independently of the value of S. 

A second class of N spins S density matrices we con- 
sider is a generalization of the Werner density matrices 
H, given by §) with 



>= 



1 5 
Af ^ 



a k \k, k, 



,k>, M= 



(16) 



k=-S 



k=-S 



The standard spin-S" Werner density matrix, corresponds 
to a_5 = ■ • • = as and it was shown [ B8J that for this 
case the exact value where quantum entanglement takes 
place is x c — 1/(1 + D 1 ^^ 1 )- For this particular case 
it was also shown recently [24|] that the entropic bound 
( |l3| ) also coincides with the above exact value xf = x c . 
The application of the Peres criterion by transposing the 
last spin give us the eigenvalues ±\aja,k\/ftf (I ^ k) or 
|«fc | 2 /A/" for the matrix a in (||) and consequently the 
bound (0) is given by 



< = (l + U JV Max{|ofo fc |}/ 



s 

E 

k=-S 



Ofc 



(17) 



where we denote by Max{|a ; *afc|} the maximum value of 
the product |a*afc| (I ^ k) of the coefficients forming the 
general density matrix (||) and ([ll]). On the other hand 
the largest eigenvalue v = Max{|ofe| 2 } / 'J\f for the matrix 
p in (pd|), give us, from (|l^) the entropic bound 



i+ ^Max {Kflfc | } 

£L-sKI 2 

a _ Et-5KI 2 -Max{|a fc | 2 } 
(D-l)M&x{\a*a k \} 



(18) 



For arbitrary values of {a k } we have a < 1 and from (jTi 
and (18) we see that the Peres criterion is in general more 



restrictive than the entropic one, i. e., x p < xf. For the 



standard Werner density matrix (a_s 



1) we 
that 
Like 

5 



■=a s 

have a = 1 and the result x^ = x~ = 1/ (1 + D N ~ 
in this case coincides with the exact value of x c p 
in our first application ( [14] ) and (15), whenever x' c = x 
these values turn out to be the exact one. Equation (|17| ) 
and (|l8|) tell us that, whenever at least one of the com- 
ponents aj, = or at least a pair of coefficients exists 
such that a k ay in (|l6|), a > 1 and these bounds are 
distinct. The coincidence of bounds x^ = happens 
only for the special case a~s = . ■ ■ = as- 

It is important to notice that interpreting (Q) with ( |l6| ) 
as the coupling of an entangled state \cf> > with the white 
spectrum environment, in order that the system stays 
entangled, any component of the environment should 
not exceed a fraction r c = (1 — x c )/(l + (D N — l)x c ) 
of the component of the engineered entangled state </>. 



3 



Using the more restrictive bound x^ ( |l7| ) we obtain 
r c = (Efcl a fc| 2 " Max{|a;a fc |})/Max{|a*a fc |}, that im- 
plies that r c decreases with the number of components 
where au 7^ in ( |l6| ) and it is more difficult to keep the 
system entangled. 

We have also studied other density matrices where the 
special entangled state \<j> > are invariant under spatial 
translation on the spin ordering as well charge conju- 
gation. In this case the results have to be carried nu- 



merically 30 1 and our results indicate concidences of the 



bounds x 1 ^ and xf as the number of subsystems N — ► 00. 

In conclusion we have shown that, although the Peres 
transposition method gives us more restrictive bounds 
than those coming from the nonextensive entropy, the 
application of both methods, that in general, at least nu- 
merically, is not a difficult problem, have the advantage 
of possibly providing, when the bounds xf and xf coin- 
cide, exact results, that are always of difficult derivation. 
Also we should remark that our analysis provides inter- 
esting results that should be useful for quantum com- 
puting, teleportation or cryptography. If the pure engi- 
neered entangled state \<j> > is the desired information 
we are processing (x = 1 in (Q)), while time runs, the 
coupling with the environment (represented by x in (Q)) 
grows (x decreases) and the decoerence effect increases. 
As long as the system stays entangled, it might be pos- 
sible to recover the desired information with appropriate 
correcting, distillation-like, procedures. Then, once we 
know the time dependence of x(t), the bounds xf and xf 
will give us the time scale where the information should 
be processed. 
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